Construction of topological defect networks with complex scalar fields 
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This work deals with the construction of networks of topological defects in models described 
by a single complex scalar field. We take advantage of the deformation procedure recently used 
to describe kinklike defects in order to build networks of topological defects, which appear from 
complex field models with potentials that engender a finite number of isolated minima, both in the 
case where the minima present discrete symmetry, and in the non symmetric case. We show that 
the presence of symmetry guide us to the construction of regular networks, while the non symmetric 
case gives rise to irregular networks which spread throughout the complex field space. We also 
discuss bifurcation, a phenomenon that appear in the non symmetric case, but is washed out by the 
deformation procedure used in the present work. 

PACS numbers: ll.lO.Lm, 11.27.+d 



I. INTRODUCTION 



\ Defect structures have appeared in high energy physics almost fifty years ago, with some of the pioneer results 

■ collected in Refs. El, Si II. Along the years, the subject has grown in importance, accompanied with an increasing 

number of investigations on kinks in one spatial dimension, vortices in two dimensions and monopoles in three 
. dimensions, among other topological defects - see, e.g. Q for an extensive discussion of some of the most important 
f*) i' results in the area. 

^ , The classical solutions which represent the defect structures can be of topological or non topological nature, and 
here we will deal with perhaps the simplest topological structures, which appear in models of scalar fields. To be 
\ specific, we will consider models of the Wess-Zumino type, described by a single complex scalar field in the presence 
' of discrete symmetry and in the more general case which engenders no specific symmetry. Some of the models 
] have been studied before in 0] - see also Q for related issues - with particular attention to the presence and 
. stability of kinklike defects and junctions, and in where the kink orbits are written in terms of real algebraic 
curves and the equations of motion are shown to be fully expressed in terms of first order differential equations of the 
Bogomol'nyi-Prasad-Sommerfield (BPS) type Q. 
• ■ The kinklike structures have been used in many different contexts, in (1, 1) and in higher space-time dimensions, 
Jli^ in particular in the form of junctions and networks of defects 0, [13, [O, f'^ 1) dimensions they are usually 
named domain walls, which can find applications in several distinct scenarios, in particular as seeds for the formation 
^ of structures in the early Universe. In this context, although the standard scenario seems to show that the presence 
• • I of domain walls has little to contribute to the cosmic evolution, it has been suggested that domain walls may perhaps 
be used as a source for the dark energy necessary to feed the current cosmic acceleratioi i IIJ . 

Another line of research has recently appeared in gravity in higher dimensions [13, Il5||. with the hope to solve 
5^ ' the hierarchy and other problems in high energy physics. In (4, 1) dimensions, the braneworld model with warped 
geometry involving a single extra dimension of infinite extent suggested in [ist has strongly impacted the subject. In 
this braneworld scenario, the inclusion of scalar fields may contribute to smoothen the brane |16|. to give rise to a 
diversity of situations of current interest, as one can see, for instance, in the recent investigations |17|. 

The present study is a continuation of a former wor k llSl. Here we will focus mainly on the deformation procedure 
introduced in [l^ , and extended to other scenarios in [20l.l2ll| . Those investigations have led us to find a peculiar and 
very interesting feature of the deformation procedure there implemented. The issue is that it is sometimes possible to 
deform a given model described by a potential containing some minima, to get to another model, with the potential 
giving rise to a different set of minima, which may increase periodically. This feature strongly suggests the possibility 
of using the deformation procedure to build lattices of minima in the two-dimensional field space. 

An interesting property of the deformation procedure is that it also constructs the kinklike solution of the deformed 
model in terms of the kink solution of the original model. Thus, in the lattice of minima we can then nest a network 
of defects very naturally, that is, as internal feature of the deformation itself. This is the idea underlying this 
paper, in which we apply the deformation procedure to investigate the generation of networks of kinklike defects for 
the deformed models, which are expanded networks. Although it is possible to start with the more general case, 
considering models with an arbitrary set of minima, we shall firstly deal with the case involving N minima in a Ziq 
symmetric arrangement. We shall consider the symmetric iV = 2, = 3 and A^ = 4 cases explicitly, and later we relax 
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the constraint to deal with three and four minima in the non symmetric case. The main reason for this is that we want 
to keep the motivation set forward in our former work [l^, where we have investigated the construction of regular 
networks. Moreover, as a pedagogical concern we believe that this route makes the problem easier to understand. 

The idea of constructing networks of defects is not new, but the novelty here relies on the use of the deformation 
procedure as a simple and natural way to generate networks. The mechanism is powerful and suggestive, and fully 
motivates the present work. To make it short, direct, we have decided to consider models of the Wess-Zumino type, 
driven by a single complex scalar field. These models are popular, of great importance and easy to manipulate, and 
so they very much help us to highlight the idea to be explored below. 

We start the investigation in Section [TTl where wc introduce the symmetric models and perform the deformation 
procedure on general grounds. In Section IIIII we illustrate the procedure with some applications, considering two 
important cases, which engender three and four minima, forming an equilateral triangle and a square, respectively. 
There we show how the deformed models tile the plane replicating the sets of minima in the entire field plane. 
We then consider other possibilities in Section IIVI and there we deal with more general models which three and 
four minima, engendering no symmetry anymore. We use the deformation procedure to get to many distinct and 
interesting patterns. The more general case allows for a new phenomenon, bifurcation, and so in Section |V] we deal 
with bifurcation, which concerns the possibility of the system to allow for two or more distinct connections between 
two given minima. This is related to the marginal stability curve, and has to do with the energy balance involving 
distinct orbits in field space, as already investigated in 0. We end the paper in Section IVll where we introduce some 
comments and conclusions. 



II. DEFORMATION OF WESS-ZUMINO MODELS 



In this Section we start with a brief summary of the bosonic sector of the standard Wess-Zumino model engendering 
the Dn symmetry. We then propose a simple but very interesting way to deform the model, to generate an infinite 
family of new Wess-Zumino like models with their defect solutions. 



A. The general case 



Let xi^T^) = Xii^^'t) + *X2(a^,i) be a complex scalar field, written in terms of the two real partners and 
X2ix,t) in (1, 1) spacetime dimensions. The dynamics of the bosonic sector of Wess-Zumino models is governed by 
the Lagrange density 

c = ld^xd''x~v{x,x) (1) 

where the bar stands for complex conjugation. We refer to these systems as Wess-Zumino, or Landau-Ginzburg, models 
if the potential energy density is determined from an holomorphic superpotential such that the potential energy 

density of the scalar field theory reads 

1, 



Vix,x) = ^W'{x)W'{x) (2) 

The interest of these models lies in the fact that all of them admit a supersymmetric version with N = 2 extended 
supersymmetry. It is also important to stress that there is a C/(l) ambiguity in the election of the superpotential: 
Wa{x) = e~*"M^(x), where e'" S U{1), produces the same dynamics. 

We shall firstly consider polynomials of degree -|- 1 in x with real coefficients as superpotentials, in this case 
W{x) = W^(x) is the same function of the conjugate complex field. The vacua manifold, the set of zeros of V{x,x)^ 
is given by the critical points of the superpotential, the N roots of the polynomial W'{x) 

^i3) (^x) = v^^\ yy'(„O))^0, i = 1,2, •••,7V (3) 
The BPS static kinks satisfy the system of first-order ordinary differential equations 

^ = % - (4) 

We use these expressions to show that the field profiles and orbits should obey 

dx dx 



dx==^^—f-, W^ix)dx-W^,ix)dx = (5) 
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Equation ([5]) (right) means that d{Wa — Wa) = such that 

ImWaixix)) = constant 
for the field orbits. These orbits arc kink orbits if they connect two vacua 

ImWa{k,) {v'^''^) = ImWaikf) {v'^^^) = constant 
This criterion sets a^''^'^ by requiring that (W{v^'''>) — e~'"' be real, or 



arctan 



Integration of (O (left) gives 



Re(M^(i;(fc))-M^(wO))) 

^ f ( dx , dx 



modTT, a^^'') = a^*''-'^ + • 



(6) 
(7) 

(8) 
(9) 



Interpretation of the meaning of this integral is reached from the identity between differential one-forms 



W:,{x)dx + W^^{x)d^=2\W\x)\ dx ^ d{Wa^ + Wo.) = 2\W'{x)\ dx 



\W'{x)\ dx = s (10) 



The kink profiles are then obtained by inverting these relations between the real part of the superpotential and the 
"length" s on the kink orbits ^ - see, e.g. Ref. Q. 
The energy of the static configurations can be written a la Bogomol'nyi in the form 



E=-ldx 



(11) 



which shows that the solutions of the first-order equations (|4]) for the kink phases a'^'^-') ([8]) have energies given by 

M(fcj)= Rel^„(.,)(w"'')) -Rel^„(.,)(w(j')) = W{v'^''^) - W{v'^^'>) 



B. Symmetric Wess-Zumino models 



The choice of the holomorphic superpotential in the form 



wix)^xi^,t) 



x,t) 



N ■ 



N en 



leads to the potential energy density 



nx,x) = :^(l-X^(x,t)) (l-x^(a;,t)) 



(12) 



(13) 



The dynamics is invariant with respect to the dihedral group = Z2X Zjy, the symmetry group of a regular polygon 
of N sides. In our case, the Z2 sub-group is generated by the transformation x ^ X ^-nd the elements of Zn are: 



X 



e « 



X , n = 1, 2, • • • , TV. Because the vacuum manifold is the set of the A^th roots of the unity 



x''''Hx,t) = v^"' = exp(27ri(fc- l)/iV), A: = 1,2,, 



the Dn symmetry is spontaneously broken to the complex conjugation Z2 sub-group at every vacuum state. 
The BPS static kinks satisfy the system of first-order ordinary differential equations 



dx 
dx 



:e-(l-x^(x)); 

These equations can also be written as 

e~'" dx e*" (ix _: 



dx 
dx 



\1-X'\x)) 



dx 



i-x 



N 



1-x^' 



(l-x'')rfx-e»(l-x'^)dx = 



(14) 



(15) 



(16) 
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Thus, the real algebraic curves which solve ((T6|) (right) 



Im 



constant 



N + 1 

are the orbits of the solutions. Kink orbits pass through two minima of the potential, so we have 



Im 



N + 1 



N 



27r 

sin(— (fc — 1) — a) = constant 



and 



Im (z;('=) ) = Im {v<^^^ ) ^ 



(N + l) 

Q,(fei) ^ -arcsin (cos {f{k + j- 2))) , k> j 
Q,(fej) ^ arcsin (cos (-^(fc + j - 2))) , fc < j 



(17) 



(18) 



(19) 



where a^^''^ — a^''^^ + tt. 

Integration of pB]) (left) gives 



X-.T0 = - 



dx 



e^^{l~X^{x)) e-»(l-x^(a;)) 



(20) 



This integral can be written in terms of the local parameter ^ = |e '"(1 — {^))\ i^i order to implicitly obtain the 
kink profiles 



Re(e— (l-x^(a;))) 



{1-X''(.m dx 



The kink energies are 



M{kj) = 



2N 

N+l 



sm{—{k-j] 



(21) 



(22) 



In [y] it was shown that this superpotential in the A/" = 2 supersymmctric Landau-Ginzburg action is an integrable 
deformation of the A/" = 2 supersymmctric minimal An series of conformal models. It was also suggested in [g the 
connection with the solitons of the affine Toda field theories - see Ref. Q for details - which can be directly 
envisaged in the above expression 



C. The deformation procedure 

We now turn attention to the deformation procedure, which will allow us to obtain new Wess-Zumino like models. 
According to Refs. [H, [1^, we will express the deformed system in terms of a new complex field (j){x, t) = 4>i{x, t) + 
i4>2{x,t), related to the original one by means of the (a priori) holomorphic function /(</>) such that 



X = /(0)=./l(<^l,02)+*/2(01,02) 



This function has to obey 



The first-order equations become 



d(j>i d(j>2 ' 



dfi 



dh 

901 



dx 



n4>) 



dx 



(23) 



(24) 



(25) 



that we choose to understand as determining the absolute energy minima associated to the "deformed" Lagrange 
density 



(26) 
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The dynamics governed by C and Cd are different, but we can define V((^, 4>) and W{4>) by 



such that the "deformed" first-order equations (PS)) are 



^=e'"W'(0); ^=e-'"W'(<^) (28) 

The BPS kink solutions for this system are obtained from the solutions of pS]) by simply taking the inverse of the 
deformation function: <f)^ [x) = /~^(x^(a;)). Thus, we can make the following relation between the deformed and 
original equations: if {^) is a kinklike solution of the original model, we have that 

lmW{x^{x)) = constant; RcW{x^{x)) = s (29) 

and so we get that (p^^ {x) — f~^{x^i^)) is kinklike solution of the deformed model, obeying 

ImW(/-i(x'^(x))) = constant; ReW{r\x'' i^))) = ^ (30) 

where a is defined by 

a = / \W'{r\x''{x)))\''dx (31) 



Alternatively, one could understand ((25)) as the first-order equations of the original model written in the form 

1 



^ = ^.f'Wn<t>)d^<f>d''^ - V{f{^), f{4>)) (32) 

This interpretation means that the original system in the new variables appears as a nonlinear sigma model with 
target space a non-compact Riemannian manifold with metric 



G^40,0)=O = G^^(0,0); G^^(,/,,0) = /'(0)/'(0)^G^^(0,0) (33) 

The merit of our approach is that we infer the kink solutions of one complicated but interesting field theoretical 
model from the well-known kinks of the associated simple system. The other point of view, in which one deforms the 
metric rather than the potential energy density, is sometimes also interesting. Despite dealing with the same model, 
the use of appropriate coordinates in field space may lead to separation of variables in the first-order equations, 
sometimes reducing its integration to quadratures; see e.g., [2^ . 

Inspired by former investigations on the deformation procedure, we select the deformation function as being equal 
to the new superpotential, that is, we choose /(0) = VV((/)). This choice constrains the function /(</>) to obey the 
equation 



/'(0)/'(0) = V^2F(/(0),/(0)) (34) 

A function / satisfying this condition assures the relation (p7)) to be fulfilled and presents the advantage of providing 
a potential for the new model which is well defined (finite) at the critical points of /(</)), i.e. the zeros of f {<!>)■ As a 
bonus, the procedure leads to a very simple expression for the deformed superpotential. 



III. DEFORMATION OF SYMMETRIC WESS-ZUMINO MODELS 



To clarify the general considerations, let us now illustrate the above results with explicit examples. We will consider 
the cases N = 2, N = 3, and N = 4. The case = 2 is simpler, and it is very similar to the deformation used in the 
first work in to get to the sine-Gordon model. The cases A = 3 and A^ = 4 are harder. The deformation procedure 
leads to the formation of junctions of kink orbits from the original Wess-Zumino kinks. Since the original non deformed 
models engender sets of minima which depict equilateral triangles and squares, respectively, the deformation will then 
naturally tile the plane, with networks of defect orbits which we name expanded kink networks. 

We will solve ^ for the Wess-Zumino model with solutions of f'{(j))^ = (-1)^(1 - (</>)), and for A = 3, 4 these 
solutions are meromorphic functions. The issue here is that the deformation function /((/)) fails to be holomorphic in 
a discrete (infinite) set of points, T, and this induces the potential energy density to acquire a countably infinite set of 
poles (the metric in the target space in the second approach above acquires a countably infinite set of zeros) . Analogous 
physical systems are described by the elliptic Calogero-Moser models (the elliptic tops in the second framework); see 
e.g., Ref. [2^. The loss of holomorphicity can be avoided by restricting the new field to take values away from the set 
r, the lattice of poles of f {(/))■ We shall then take C/F as the 0-field space. 



A. The case N = 2 

111 this case, we define the field x as a function of the new field <f) in the form x — f{<P)- With this, we can use 
to rewrite 

wix)^x-lx'; nx,x) = ^(i-x')(i-F) (35) 



as 



W{f) = /(^) - Ifi^y, Vif) = i(l - - Pm (36) 



The deformation function f{4>) is the new superpotential if we impose 



/'(0)/'(0) = ^Jii-pma-pm (37) 

The particular choice /(0) = sin((^) complies with ([57]) and leads to the deformed system defined by 

W(0) ^ cos{(f>) ; V(</.,0) = cos(<?!)) cos(^) (38) 
which is the complex sine-Gordon model. Here the first-order equations are 

^ = e*" cos(^); ^ = e-'" cos((/)(x)) (39) 
ax ax 

The solutions for a = 0, tt arc given by 

[x) = ±gd(x) -I- 2n7r = ± arcsin(tanh(a:)) + 2n7r (40) 

where gd stands for the Gudcrmannian function [2^, and n is an integer. Here the kinks are analytic solutions and 
the superpotential is holomorphic. 

B. The case iV = 3 

In the = 3 case, we have 

W{x)^X~\x\ V{x,x) = l{i-x'){l-t) (41) 



and putting x = /(0)i we rewrite this formula in the form 



As stated in (|34p . the deformation function must then satisfy 



/'(0)/'(0) = ^J{i-pimi-pm (43) 

We choose in particular the holomorphic solution of (|43p which satisfies the separated equations 



f'W' = fW' - 1 ; /'(0)2 = fi^r - 1 (44) 

The solution of (j44|) . henceforth a solution of (|43|) . is the equianharmonic case of the Weierstrass T' function 

>V(<?!)) = /(^) = 4^-p(4-^(/);0,l) (45) 

Recall that the Weierstrass V function - see [l^l - is defined as the solution of the ODE 

iV'{z)f ^AV''{z)-g2riz)-g3 (46) 
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The Weierstrass Viz; g2, 93) elliptic function and its derivative that solve the differential equation above are doubly 
periodic functions defined as the series 



m.n 



^'W = ^-2V- ^— (47b) 



with m,n G and + 7^ 0. Therefore, the deformation function is, up to a factor, the Weierstrass V function 
with invariants 32 = and — 1, and we denote it by Voiiz). This function is meromorphic, with an infinite 
number of poles congruent to the irreducible pole of order two in the fundamental period parallelogram (FPP). Thus, 
we suppose that the 0-field takes values away from the set of points Fa in order to make the new superpotential 
holomorphic in the = 3 case. 

A brief reminder of the essential properties of Voi and Vqi is the following: 

1. 7^01(4"^^) is a single-valued doubly periodic function with primitive periods: 2loi and 2uj3. Defining uj2 = 
43r'^(l/3)/47r, the primitive half-periods are 

-i=-2(^i-*^^; u;s = u;2(^l+^^^ (48) 

These periods determine the FPP. Note that only two of the half periods are irreducible: lu2 = coi + UI3. 

2. 7^01 (4~^0) has only one pole of order two at = in the FPP. 

3. The values of 7^01(4" 3?!') at the half-periods are 



7'oi(4-^wi) =4-3; 7?oi(4-3w2)=-4-3 --hi— ; 7?oi(4-3a;3) = -4-3 j (49) 





4. The zeros of the derivative 7^01(4 ^4') iii the FPP are at the half-periods of Vqi 

r^^{A-iuji) = 7'^i(4-3^2) = 7'^i(4-^^^3) = (50) 
and the origin is a third-order pole of 'Pqi{z). 





Figure 1: (Color online) The symmetric case A'' = 3. 3D graphics of the potential V{(f>,(l>) (left panel), and of —V{4),(j>) near a 
point in the lattice Fa (right panel). Note that in the right panel the zeros are now maxima. 



With these ingredients we write the deformed potential 

1 



V(0,0) = ^^J{^~4V|,{A-^)){l-iV|,i4-^)) = -7'^,(4-3</,)Poi'(4"3, 



(51) 



which is depicted in Fig. 1. 
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The new potential is doubly periodic with an structure inherited from the "half-periods" of V. The set of zeros of 
the potential in the FPP (see Fig. 2) has three elements 



1 .\/3^ 



0(3) ^ = 43 



The set of all the zeros of V form a lattice (see Fig. 3) which tile the entire configuration plane 

1. 



(m,n) \ 9 



+ - + \/3 z(m — 71 + - 



(m,n) 



CJ2 



^TO + 71 + 1 + A/3i(TO — n)^ 

The values of the superpotential at the minima are 



(52) 

(53a) 

(53b) 
(53c) 

(54) 





Figure 2; (Color online) The symmetric case = 3. 
connecting the zeros of the potential (right panel). 



Zeros (red) of V and points (yellow) of the set Fa, and kink orbits (blue) 



In sum, the potential obtained from the deformation procedure has the same zeros in the FPP as the original 
model. Besides, one pole arises at the origin due to the meromorphic structure of Pqj; see Fig. 2 and 3. However, 
this structure is infinitely repeated in the deformed model, according to the two periods uji and determining the 
modular parameter r = ws/wi = — 1/2 + i\/3/2 of the Riemann Surface of genus 1 associated with this T'-Wcierstrass 
function. As an aside, we note that the modular parameter r = 1/2 + i^/2>/2 gives the same Riemann surface because 
r = [ar + h) / [ct + d) where 





















is an element of the modular group SL{2, Z). 

Contrarily to the deformation function chosen in the case = 2, for = 3 the Weierstrass V function is not an 
entire function, that is, it is not holomorphic in the whole complex plane C. F = flmn = 2muji + 2nL02 {m, n E Z) is 
the lattice of points of V and V which are accordingly meromorphic functions. Thus, we suppose that the new field 
take values in the space C/F3 to avoid the loss of holomorphicity. This point of view is very close to consider the 
genus 1 Riemann surface C/~F of modulus t = 103/0^2 minus the origin (the FPP with the edges identified pairwise 
minus the origin) as the 0- field space. Keeping, however, the infinite copies of this space contained in C/F3 gives a 
richer kink structure. 
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B.l. Network of P-kink orbits 



We shall compare the 7^-kink orbits with the orbits of the original = 3 polynomial Wess-Zumino model. If {x) 
is a iV = 3 solution of and (HO]) then ^'^{x) = A^Vq^{\^^x^ {^)) solves 



where 



Ime^*"4*7'oi(4"30-f^(a;)) = constant ; Re e"*" 4* 7'oi(4"*(/)" (x)) = cr 



dx 



(55) 
(56) 



-26^2 —0J2 



ii}2 2(^2 36>2 

• • • • 



( 


» • 










-2co2 /\ 






1 V ^ 



Figure 3: (Color online) The symmetric case = 3. Lattice of zeros (red) of V{(j), cj>) and points (yellow) of the set Fs, and the 
network of kink orbits (blue) in the lattice of minima (right panel). 



Because of the relations 
the same values of a as in the non deformed case, 



a 



arctan 



sin^(fc-l)-sin^(j-l) 
cos2f(fc-l)-cos2^(j-l) 



(57) 



(58) 



give the deformed kink orbits. 

There are three types, which we show below. 



Type (13), non deformed 

The condition ImWaix^^'') — IinWQ.(x^^'') is satisfied only for a''^^' = tt/6 (for antikinks) and a^^'^' = 7tt/6 (for 
kinks). What is called kink and what is antikink is a matter of convention. Our convention is that kink/antikink 
orbits run clock/anticlockwise in the (W, W) plane. The orbits obey 

_^<Reiy,(^^)<^ ; ImM^.(x^) = -^ with ReW^.(x«) = ^; ReW^ix^'^) = 

and 

-^<ReTy^(x^)<^; W-(X^) = ^ with Reiy^(xW) = -^; ReW^^ (x^^)) = ^ 

Type (13), deformed 

As mentioned above, the condition ImyVQ,((/)|j2 „)) = ImWa (^^^j^, ^,^) is again satisfied only for 

^ 7r/6 and 

a(i3) ^ 77r/6. The orbits obey 

- — < RcW^U^) < — ; ImWiLU^') = -- with ReWiL(</)p^ J -— ; ReWiLU\^\ ,J = — 
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and 

-^<RcW-(0^)<^; ImW^(^^0 = i with RcW^(0W „,j) = -^; RcW^(0S,„p = ^ 
The nearest neighbor type (1) and type (3) minima are connected by the orbit following the sequence 

^(1) ^ ^(3) ^ 0^(1) ^ ^ ,(1) 

'{m,n) ^{m,n) ^ (m-{-l,7i) ^(m+l,n) ^(?7i+2,n) 

See Fig. 3. 

The other two cases (23) and (12) follow similarly. Here we just add that for the case (23), (m, n) and (m', n') are 
restricted to link nearest neighbor type (3) and type (2) minima along the orbit. The sequence is 

^(2) J,(3) w!,(2) J,(3) i{2) 

^(m,n) ^(m,n) ^{m,n+l} ^(m,n+l) ^(m,n+2j 

Also, for the case (12) we have that (m,n) and {m',n') must be chosen according to the following sequence 

^^(2) ,(1) ,(2) ,(1) 

in order to connect nearest neighbor minima of type (1) and (2). 

We end the case N = 3 collecting the corresponding energies. We have that the defect energies for the original 
non deformed Wess-Zumino model are given by M = 3^3/4, for kinks and anti-kinks for all the three sectors, with 
(f^j) = (12), (23), and (13). For the deformed model, the energies of the T'-defccts are given by = a/3, for the 
same cases. 

C. The case N = A 

In the = 4 case we deal with 

W{x) = X-lx\ Vix,x)^l{l-x'){l~t) (59) 

Thus, putting x = /(0)i we have 



w{f)^fi4>)~-fH^), V ^ -{1 - fma - f'W) (60) 



The special deformation function must satisfy 



Arguing like in (j43p . we choose the holomorphic solution that satisfies the separated equations 



/'(0)^ = 1-/(0)% /'(0)2 = 1 - /(0)4 (62) 

The solution of (|62[1 . henceforth a solution of ([M]) . is the elliptic sine of parameter fc^ = —1, the Gauss's sinus 
lemniscaticus 

W(0) = /(0) = sn(0,-l) (63) 

As the derivative of the Jacobi elliptic sine is sum' = cnii dnii, the identities 

cn2(0,-l) = l-sn2(0,-l); dn2(0, -1) = 1 + sn2(0, -1) (64) 

show the solution of (pTjl very directly. 
The deformed potential reads 

V(0, 0) = iv^(l-sn4(0,-l))(l-sn%0,-l)) = i|cn(0, -l)P|dn(0, -l)p (65) 

which is depicted in Fig. 4. The superpotential is given by Wq(0) — e^*" sn(0, —1) 
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Figure 4: (Color online) The symmetric case A'' = 4. 3D graphics of the potential V((/>, (fi) (left panel) and of —V{4>, (p) near four 
points of the set r4 (right panel). Note that in the right panel the zeros are now maxima. 



The new potential is doubly periodic with its structure inherited from the "quarter-periods" K{~1) = uJi/A and 
iK{2) = u)2l^ of the twelve Jacobi elliptic functions; see Fig. 5. Here K{—1) w 1.31103 is the complete elliptic 
integral of the first type, a quarter of the length of the lemniscate curve in field space: {(t)\ + 4>2)'^ = 4)\ ~ (t>\. 
K{2) « 1.31103 — il. 31103 is the complementary complete elliptic integral of K{—\). 

The set of zeros of the potential in the FPP are (j)'^^^ = [Ji/4, (Z)*^^^ = /4, (/)('^) = — a;i/4, (/)^"') = — ia;i/4, whereas 
the set of all the zeros of V form a quadrangular lattice in the whole configuration space. They are given by, explicitly 

4>M = T (^^^^^ + «) + ! + ''an) ; 4>fl^) = ^ (2(2m + n)+ z(2n + 1)) (66a) 



4 

(m,n) 4 ^ ^ ' ' ^(m.ri) ^ 



- '^{2{2m + n)-\ + i2n)- 41 . ^ [2[2m + n) + i[2n - \)) (66b) 



because 



Thus, 



cn(</'[^,„).-l)-dn(0£,„),-l)=O, fc = l,2,3,4 (67) 



Wc(<^(i') = e-'", Wa(0(2)) = ie-'", Wa(0(3)) = -e-*", W„(0('^)) = -ie"*" (68) 

since sn[±/^(— 1), —1] = ±1 and sn[±iX(— 1), —1] = ±i. 

The modular parameter of the associated genus 1 Riemann surface is t = i/\ [2]/X[— 1] = 1 + i. Identical Riemann 
surface is associated to the lemniscatic case, 92 = 1, .93 = 0, of the Weierstrass V function. Like in the former case, 
however, the Jacobi elliptic sine is not an entire function, and so we restrict the new field to live in C/r4 (where r4 
is the set of poles of /(</>) in this case) in order to make the superpotential holomorphic. 



C.l. Network of sn-kink orbits 



We shall compare the sn-kink orbits with the orbits of the original = 4 polynomial Wess-Zumino model. If }^{x) 
is a solution of (|^ then (f>^ {x) = sn^^(x^ (-i;), — 1) solves 

Ime-'"sn(</)^(a;),-l) = constant; Re e~*"sn(0-f^(x-), -1) = cr (69) 

where 

a = J |cn[(/)-^(a;),-l] •dn[</.^(a;),-l]|^ (70) 
The same values of a as in the non deformed case give the kink orbits 
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selects 

■sinf(A:-l)-sinf(j-l) 



q;('^-') = arctan 



cosf(A:-l)-cosfO'-l) 



modTT (72) 



as the angles for both the original and deformed kink orbits. 
There are four types, which we show below. 



Type (12)/(34), non deformed 
In this case we have lmWc{x'^^'^)=lmWaix'-^^) and ImM^„(x(3))=imVl^^(^(4)) Q^jy foj. 

a = 37r/4 (kinks) or a = 77r/4 

(antikinks). The kink orbits obey 

<ReW3^(x^)<^; W-(x'^) = -^ with or 
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whereas for the antikinks 



ReWs^ (;^(3) ) = 2^ ; ReW^ (^H)) = 



o/o o/o 9/0 Rel^-(x(^)) = ¥; Rciyif(x(^)) = -¥ 

_fXf <ReT^^(;^^)<f^; linW/^(x^) = i^ with or 

^ ' ^ ' ^ ReM/i.(x(3)) = -¥ ; RcM/2f(x^^^) = ¥ 

Tj/pe (12)/(34), deformed 

Here wc have ImW„(0j^ „)) = ImW,(0[^, _„,)) and ImW,(0[^| „)) = ImWo(<^[;2',«')) ^^^^ " = ^^/^ ('^™^') 
a = 7tt/4 (antikinks). The kink/antikink orbits obey 

/K /K fK ReWa^tiJ^^ ) = -^- Re>V3x(0p^ ,0 = ^ 

/o /o ./o — vr(m,n)^ 2 ' — VT'(m'^„')y 2 

- — < ReW3. (6^) < — ; ImWs. (0-^) = -— with or 
and 



ReW^(<^[^, „,)) = ^ ; ReW^(0[:^,„)) = 



^2 ^/2 ^/2 R'^Wif Cv„,) - 

- — < ReW7^((/)-^) < — ; lniW2j^{4)^) = — with or 

ReWz.(^g,,„,) = -^; Re>Vz.(C„)) = ^ 

where {m,n) and {m',n') are restricted to link nearest neighbor type (1) with type (2) and type (3) with type (4) 
minima along the orbit. The sequences are 

<f(m,n) ^ '?(m,n) ^V'(m-1,„+1) ^ <?'(m-l.n+l) ^ 'f'(m-2,n+2) 



and 



6(3) ^^(4) ^J,(3) ^ .(4) ,(3) 



See Fig. 6. 

The other three cases, (13), (14)/(23), and (24) follow similarly. Here we just add that in the case (13) we have 
that (m,n) and {m',n') are restricted to link nearest neighbor type (1) and type (3) minima along the orbit, and so 
they must be chosen according to the following sequence 

In the case (14)/(23) we have that {m,n) and {m',n') are restricted to link nearest neighbor type (1) and type (2) 
minima respectively with type (4) and type (3) minima along the orbit. Therefore, the sequences are 

'^(m,n) '(m,n) '{7n,7i-\'l) ^(rn,n+l) ^ {771 ,n-\'2) 
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1 iKi 




-K 

1 -iK' 


f 




Figure 5: (Color online) The symmetric case N — 4. Zeros (red) of V and points (yellow) of the set r4, and the two (blue and 
black) possible orbits connecting the zeros of the potential. 



• 'A* 




Figure 6: (Color online) The symmetric case TV = 4. Lattice of zeros (red) of V(0, (j>) and points (yellow) of the set r4, and the 
networks of kink orbits (blue and black) in the lattice (right panel). 



and 



/,(2) ^ A3) ^ A2) ^ A3) ^ ,(1) 

(m,n) '{ni.n) ^(m,7i+l) ^(7n,n+l) ^(?n,n+2) 



In the case (24) we have that (m, n) and (m', n') are restricted to link nearest neighbor type (2) and type (4) minima 
along the orbit. Thus, they must be chosen according to the foUowing sequence 

J^) ^ ^(2) , (4) , (2) 

We end the = 4 case collecting the energies of the defect structures. For the non deformed kinks we get 



M(13) = i\/(24) = - ; A/(12) = Af(34) = — ; M(14) = M(23) = — 
5 5 5 

whereas the energies of the deformed sn-kinks read 

X(13) =X(24) = 2; M{12) = M{M) = V2; 7W(14) = 7W(23) = \/2 
with M(fcj) = M(jk) and M{ki) = M{jk). 



(73) 



(74) 



IV. DEFORMATION OF ABRAHAM-TOWNSEND MODELS 



Let us now move on to the case where the original model engenders no specific symmetry. This study is inspired on 
Ref. 0, in which Abraham and Townsend consider some interesting situations, guided by more general superpotentials, 
which develop no specific symmetry. Similar potentials were also considered in but there the investigation was 
mainly on the symmetric case. In [7|, however, the focus was on the non symmetric c ase, as a basic model underlying 
the study of intersecting extended objects in supcrsymmetric field theories - see also [25, which deals with the forces 
between soliton states in the same model. 
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The absence of symmetry makes the investigation harder to follow, but it is still of current interest since it leads to 
more general possibilities, bringing some new effects into the game and allowing for the presence of irregular network 
of defects. We follow as in the former section, and below we consider models which develop three and four minima 
with no specific symmetry anymore, using T3 and r4 as the set of poles of /(</>) as before, but now in the asymmetric 
cases, with iV = 3 and = 4, respectively. 



A. Irregular network of P-kink orbits 



Let us start with the simplest case, in which one considers a class of models that engenders three minima. Here we 
deal with the superpotential 



(75) 



where S = Si +162 is a complex coupling constant which parametrizes the family of models. In this case, the potential 
is given by 



v{x)^l\s-x\'\i-x'\' 



(76) 



This choice leads to the following set of minima: two real minima which are fixed to be at = — 1 and V2 — 1, and 
a complex minimum at ~ 5, which may move in the complex plane for different choices of the complex parameter 
S. This is the most general case with three arbitrary minima in the complex plane, since one can always choose 
the straight line joining two vacua as the abscissa axis, crossing the perpendicular ordinate axis through the middle 
point between the vacua, setting the distance between them to be 2 by an scale transformation. The values that the 
superpotential Wa{x) ~ e^'"VK(x) takes now at the vacua are 



From these expressions we obtain 

W'a(l)-M^a(-1) 



■de- 



W^iTl) - W^{6) = ^(5 t3)(5± l)3e 



(77) 



(78) 



Therefore, the angles of the kink orbits are (mod tt) 



v(12) 



arctan 



Rc6 



v(13) - 



= arctan 



Im(((5-3)((5 + l)3) 



Re((,5-3)((5+l)3) 



v(23) = 



= arctan 



lm{{5 + 3)i5-lf) 
.Re((<5 + 3)(<5-l)3) 



and now the kink energies are given by 



Af fl2) = -f^^ + 52)1/2 M{23) = ^ (((<5i + 3)2 + SMSi - I? + Sir)' 
3^' M(31) = ^(((<5r 3)2 + <5i)((<5? + 1)2 + ^1)3)^ 



1/2 

/2 



(79) 



with M(fcj) = M{jk). Note that the three masses M(12) = M(13) = M(23) = (4/3)V3 for S = ±iV3, a value of 
the parameter for which the vacua lie at the vertices of an equilateral triangle, leading us back to the symmetric case 
which engenders the Z3 symmetry. 

We now go to the deformation procedure, changing x ~* f{(l>)i with /(</>) = VV'((/)). In the present case, we get that 
f{(j)) should obey 

and this now gives, in the specially simple case which we have already considered in the former section, 

f{4>? = 5- f -5f + f (81) 
This is again the Weierstrass equation, and if we define z = 4~3 and / = 4^'i/' + (5/3we find 

^'{zf = - g2^P - 33 (82) 
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where g2 = 43 (1 + and = -(2(5/3)(l - S'^/9). The solution for /(0) is then given by 



/(0)--5 + 437'(4-3^;52(5),53('5)) 



(83) 



The half-periods loi and lu^ of the Weierstrass V function are obtained from the invariants g2 and 53 by means of the 
equation 



52 = 60 ^ ; 33 = 140 ^ 1) 

(m,n) (m,n) 



6 

mn 



(84) 



where = 2mu!3 + 2nuji, with m,n G Z. We also have A = ^ 27.g| = 4((5^ — 1)^ and ei = —43(1 + (5/3), 

62 = 43(1 — (5/3), and 63 = 432(5/3 are respectively the discriminant and roots of the cubic equation which appears 
from the right hand side of ((82|) . 

It is interesting to note that the modular parameter t{5) — lO'^{5) / loi{5) of the genus 1 Ricmann surface associated 
to this P-function depends on 5. Thus, variations of 5 correspond to motions in the Ricmann surface moduli space. 

The deformed model is governed by the potential 



V((/)) = ^|7''(4-^;,g2(<5)),ff3((5)p 



(85) 



It has zeros in the FPP at V'{loi) = V'{ijJ^) =■ V'{loi + lu^) = 0. Thus, the vacua of the deformed model arc the 
constant field configurations 



=4*(a;i+an„); (/)g „) =43(tJi +W3 + f]^„); = 4^ (^3 + f^mn) 



The values of the superpotential 



^M) = - + 43 7'(4-3 0;g2((5),53('5)) e-'^ 



(86) 



(87) 



at these vacua are 



g + 2361 I 6 



W7(Cn))=(^ + 2*^^^^3(-3))e-'^ 



Therefore the angles of the deformed kink orbits are (mod tt) 

Im(5 



-+2t62j 6-*^ =6-' 
^+2363^ 6-''^ =(56- 



7(12^=0, 7(i^)=arctan 
and the kink masses become 



Rc((5 + 1) 



v(23) 



arctan 



Im(5 



Re(,5- 1) 



(89) 



X(12) = 2, M{li)^{{5i + lf + 5lY, M(2'i)={{5i-lf + 5l) 



2\ 2 



(90) 



Note that for 5 = ±i-\/3 the three masses are equal, M.{12) = A^(13) = A^(23) = 2, corresponding to a regular 
triangular lattice of minima. The sequences of minima connected by the kink orbits in these families are 



,(2) 
''^(?n-l,n+l) 

n 

(m,n— 1) 
,(3) 



V^(m-l,n + l) 
,(2) 



In Fig. 7 and 8 we plot the potential, minima and network of kink orbits, respectively, for the specific value of the 
complex parameter 5 = 1 + i. Comparison of these figures with figures 2 and 3 shows to what extent the complex 
parameter induces irregularity when it differs from 
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Figure 7: (Color online) The case of three asymmetric minima. 3D plot of the deformed potential —V{(j>, 4>) for S — 1 + i near 




Figure 8: (Color online) The case of three asymmetric minima. Plots of the zeros (red) of the potential and points (yellow) of 
the set Fs, and the networks of kink orbits (blue) connecting the zeros for 5 = 1 + i (right panel). 



B. Irregular network of sn-kink orbits 



We consider now the superpotential 

Wix)=-6ex+ li6+e)x^+ l{5e^l)x'~ J('5 + e)x'+ ^x' (92) 

where 6 and e are two complex coupling parameters which control the model. With this polynomial of fifth-order W 
we get the potential 

V{X, X) = ^|1 + xHl - X\'\e - X\'\5 ~ X? (93) 

displaying the four minima = — 1, W2 = e, 1)3 = 1, and W4 = b. We notice that two of the minima are again fixed at 
the values ±1, but this docs not restrict generality of the procedure for the same reasons as before. 
Again, we choose the deformation /((/>) = W((?!)) determined from the specially simple case 

/'(0)' = (/ + l)(/-l)(/-'^)(/-e) (94) 



This equation can be written as the elliptic sine equation y' — ■\/(l — ?/^)(l — fc^y^) if we follow the usual 
procedure [2^, in which we write the product (/^ — 1)(/ — (5)(/ — e) as the product of the two factors 
(A+(/ - a.f + A_(/ - and (B+(/ - a^f + - a^f) , where 



2l V(l-'52)(l-a / ' H 2^(i_52)(l-e2) I ± 6 + e 



(95) 
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Wc then make the homographic substitution p = {f ~ a^)/{f — «+) to obtain 

dy 



V(l-y2)(i„fc2y2) 



y(0) = sn(0, k ) 



(96) 



where y = ^—A^ /A^ p, ij.' = (a_ — 0, and fc^ = Therefore, we finaUy find that 

^ Q+ - a- y/-A_/A+ snji/j, k^) 
l-^-A_/A+ sn(^,fc2) 

is the new superpotential from which we can construct the potential of the modified model in the form 

1 {a+ -a_)2^A_S_cn('0,fc2)dn(?/', fc2) 



(V-^-M+sn(V',fc2)-l)2 



(97) 



(98) 



The above investigation is very general, and the presence of the two complex parameters S and e make the illustra- 
tions awkward. For this reason, we shall restrict ourselves to the simpler case where e = —d and only a single complex 
parameter is free. Thus, 



w{x) = s\ i(i + s')x' + \x' ; Vix.x) = x'\'\s' - x'\' 

6 2 



(99) 



are respectively the superpotential and potential of the original model. The four minima of V are at the points 
vi = — U3 = —1, V2 = — W4 = —S in the x-complex plane. 
From the values of the superpotential at the minima 



we obtain 



Wa{6)-Wa{~6) = ^<5^^(1 - ^)e-'" 
d 5 



The angles of the kink orbits are (mod tt) 

"Im(l - 55^) 



Wa.{±l)-Wa.i±S)=±- 



2 fS^^-l 



5 

2/(5^ + 1 



6\S + 1) e"'" 



a 



(13) 



a 



(24) 



arctan 



arctan 



Re(l - 5,52) 
lm{6^{S'^ -5)) 



Q,(i2) ^ ^(34) ^ arctan 



_Re((53((52-5)) 
The corresponding energies are given by 



«(23) ^ «(14) 



arctan 



Im(,52_;^_5^2(^_;L)) 

Re{S^ - 1 ~ 5S^{S - 1)) 
lm{S^ + 1-55^(6 +1)) 



Re((52 + 1-552(J + 1)) 



M(13) = -\l/5-5^\; M{12) = Af(34) = -|l/5 + (5 - l);^^ - 5^51 



Af(24) = -|((5V5- 1)(53| ; M{23) M{A1) ^ -\l/5 - [6 + 1)S^ + 6^51 



(100) 

(101a) 
(101b) 

(102a) 
(102b) 

(103a) 
(103b) 



plus M{kj) = M{jk). Note that for (5 = i, il/(13) = Af (24) = |, Af (12) A/(34) = Af(23) = Af (14) = |\/2 and we 
recover the case of a regular square. 

We now deform the model, choosing f{(j>) ~ W{(j>). As before, we consider the specially simple case, which gives 



/'(0)' = (1 + ./)(1 - f )(S + f )(S -f)=5^-il + 6'^)f{c^f + fic^Y 



(104) 



We compare this with the elliptic Jacobi sine equation to find f{(j)) — 5 sn(0,(52) as the solution of equation p04p . 
Thus, the deformed superpotential and potential are given by 



W^(0) =5sn(0,52)e-»^; V(0,(/)) = ^\S{' \cni<j>, 6-')\-'\dni<j>, S')\^ 



(105) 
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The zeros of the potential form the lattice of vacua at the constant values of the field 

1 



(m.n) 

x(3) 
(m.n) 



1 



-(Wl + UJ2) + ^ri 
4 

(wi + UJ2) + ^ mn 



where the periodicity is determined by the quarter periods of the Jacobi elliptic sine 

1 



From the values of the superpotential at the minima 



-If 



a(4) 



(2) 



we derive 



(4) 



a(2) 



The orbit angles and the energies of the deformed kinks are given by 

"Im(l - S) 



(106a) 
(106b) 

(107) 

(108) 

(109a) 
(109b) 
(109c) 



7(^3) = 0; Al(13)-2; 7(12) = ^(34) 



arctan 



.(24) 



arctan 



Red 



Re(l - 5) 

M{24) = 2\S\; 7(^3) 7(14) ^ ^^^.^^^^^ 



7W(12) = 7W(34) = \1-S\ 
'Im{l + 5) 



Re(l + 6) 



A/((41) = X(23) = \l + 5\ 

to the 



with 7^=^) = jW^^ -\- TT and M.{kj) ~ M.{jk). The vacua connected by the kink orbits are organized accordin; 
following sequences 



(1) 

{m,n) 




^ (m,n) 


<^ 






,(3) 
^(m+2,7i+l) 


(2) 

(m.n) 




-/-['^ ^ 

^ (m.n) 




,(2) 

V'(m,n+1) 




'f^(m,ri+l) 


(1) 

(m,n) 




T^(m,n) 




V^(m+l,n) 


^ 


,(2) 
'''(m+l.Ti) 


(3) 

(m,n) 




'/'['^ ) 

^ (m,n) 




,(3) 
'P{m+l,n) 


<=> 


,{4) 

''^(m+l,n) 


(1) 

(m.n) 




^ 

^ {m.n) 


<^ 


^(m+l.n+l) 


<^ 


^(m+l,n+l) 


(2) 

(m.n) 




^ {m.n) 


<^ 


,(2) 
'^Im.+l.n+l) 


<^ 


,(3) 
'?'(m+l,ri + l) 



To illustrate the investigations, in Fig. 9 and 10 we plot the potential, minima, points in the set F, and orbits of 
the kinklike configurations. 



V. BIFURCATION 

All the models which we have been studying so far present an interesting feature, which we now explore. It concerns 
the fact that they have three or four minima. Thus, if we choose two minima arbitrarily, it may be possible that they 
are connected with two or more distinct orbits. When this happens to be the case, we say that the system develops 
a bifurcation, since one can go from a given vacua to another one, following two or more distinct kink orbits. This 
possibility is directly related to the balance of kink energies providing an upper bound for the fusion of two of the 
kinks in a single kink of a third type. Such a process is energetically possible - and the outgoing kink stable - if, 
given e.g. three minima k,j, I, the energy of the (kj) kink is lower than the sum of the other two kink masses 



M{kj) < M{kl)+M{lj) 



(112) 
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Figure 9: (Color online) The case of four asymmetric minima. 3D plot of the deformed potential —V{(j>, (j)) near four points of 
the set r4 for S = 1 + i. 




Figure 10: (Color online) The case of four asymmetric minima. Plots of the zeros (red) of the potential and points (yellow) of 
the set and networks of kink orbits (blue) connecting the zeros for 5 — 1 + i (right panel). 

All the kink masses depend on the N ^2 complex parameters, indicating the arbitrary positions of — 2 vacua, since 
the other two minima are fixed at the ±1 points in the x-complex plane. A bifurcation occurs when the inequality 
becomes equality, since we can go from k to j following the direct k j path, or then visiting / through the path 
k I ~> i with the same energetic cost. In the case A^ = 3 it is shown in in the search for intersecting domain 
walls, that this possibility indeed happens - see also Ref. [1^. In the generic case A^, there is a sub-manifold of real 
dimension 2N — 5 of the parameter space characterized by the equation M{kj) = M{kl) + M{lj). We shall refer to 
this sub-manifold as the marginal stability variety because some irreducible component of it is a boundary between 
two regions of the {N — 2)-dimensional complex parameter space; one region where M{kj) < M{kl) + M{lj) (the 
(kj) kink is stable and cannot decay to (kl) and (Ij) kinks), and the other region where M{kj) > M{kl) + M{lj) (the 
(kj) kink is unstable decaying to the {kl)-{lj) kink combination). In fact, things are slightly more complicated and 
a bifurcation occurs at the marginal stability variety, with the kink orbit going to infinity beyond this point, with 
Mikj) becoming divergent, breaking the direct connectivity between the vacua k and j. 

Since the energy in the topological sector is controlled by W, we immediately see that bifurcation appears if and 
only if at least three minima gets ahgned in W space, that is, iff |M^(A:) — VF(j)| = |VF(A:) — VF(/)| -I- \W(l) — W{j)\, with 
W{1) in between W{k) and W{j). This is the general condition for bifurcation, and below we use it to investigate the 
models introduced in the former section. 



A. The case of three minima 

Here we consider the case of three minima, with the model being described by the complex parameter 5. We 
notice that bifurcation cannot appear in models which engender the Z3 symmetry, since the symmetry implies that 
M(12) = A/(23) = M(31), henceforth M(12) < M(23) M(31), etc. In the Abraham- Townsend model, the marginal 
stability curve is characterized by the alignment of the three minima in the {W^ W) plane, i.e., it is the curve in this 
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plane for which a^^^^ (S) 



These identities hold if 



Im((5 



(13) 



(5) mod TT, or, 



Re((5-3)((5 + l)3) 



Im^ 



Iin((^ + 3)(^~1)3) 



{d -S){3 + SS{S^ +P+SS-6))^0 



(113) 



(114) 



which is the algebraic equation that characterizes the marginal stability curve. We note that the curve has two 
irreducible components. 

One component is the real axis in the (5-plane: S2 = 0. For real values of the parameter it happens that all the three 
minima are aligned in the real axis of the complex x-plane, and the system cannot develop bifurcation. In the two 
very special cases with Si = ±1, we have only two minima, and so only one kink orbit remains in the system. For 
i5i < —1, there exist only two kink orbits, which we label (31) and 12, because the vacuum vi sits on the real axis in 
between V3 and U2- For —l<5i< 1, there exist the two kink orbits (13) and (32), and for 5i > 1, there exist the two 
kink orbits (12) and (23) for similar reasons. 

Things arc more interesting for the other algebraically irreducible component. The quartic curve 



3-6SI- 66l + 25l5l + ?,5\ - 5l = {5l + i - 5l + 2^ 5f - 3(5^ + 3 )((5^ + 3 - (5? - 2^ 5f - 3(5^ + 3 ) = (115) 

which allow for each one of the three possibilities, M(12) = 7\/(23) + M(31), or Af (23) = M(31) + A/(12), or yet 
Af(31) = A/(12) + M{2?>), depending on the specific value of the complex parameter (5, is the true boundary between 
regions in the 5 plane where there exist three or two kinks and the phenomenon of bifurcation takes place. This case 
is fully studied in Q, and below in Fig. 11 wc plot the curves of marginal stability in the complex 5 plane. 




Figure 11: (Color online) The curves of marginal stability in the case of three minima. The set of (red) points illustrates the 
positions of the minima and the kink orbits for distinct values of the parameter 5 in the related regions. 



The condition for alignment in the deformed model in the W plane 7 



(12) 



7 



(13) 



7(23) is 



= 



Im(5 



Im(5 



Rc((5+1) Rc((5-1) 



(116) 



There is only one component, 82 = 0, but wc have already seen that in this case there is no bifurcation anymore. 
In fact on the real axis the kink masses are A^(12) = 2, A^(13) ~ \5i + 1|, and A^(13) = \6i — 1|, such that: 
7W(13) = M{12) + M{23) if 5i > 1, M{12) = M{13) + M{23) if -1 < 61 < 1, and M{23) = 7W(12) + M{13) if 
(5i < 1. If (52 7^ the mass of one of the deformed kinks is always lighter than the sum of the masses of the other two 
kinks because of the triangle inequality, the sum of the lengths of two sides of a triangle is greater than the length 
of the third side. Thus, we notice that the deformation procedure wash out the marginal stability curve, leaving no 
room for bifurcation in the deformed models which we are dealing with in this work. In Fig. 12 we illustrate the 
N = 3 case with several distinct possibilities for the complex parameter S. 
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Figure 12: (Color online) The bifurcation curves and some related illustrations with (green) points in the first column repre- 
senting distinct values of the complex parameter S in the case of three minima. The broken (red) lines in the second column 
indicate how the orbit goes to infinity, leading to divergent contributions which should be discarded. 



B. The case of four minima 



Let us now consider the case of four minima, with the model being driven by the two complex parameters 6 and 
e. The general case is very complicate, so we move to the simpler case in which we use e = —(5, leading to a single 
complex parameter S. 

The condition for alignment of all the minima in the W plane now is 7'^^) — 7(12) _ ^(23) _ 7'^^'^), modTr (note 
that because the vacua in the pairs (vi^vs) = (1,-1) and (1^2, "^4) = (— <5, 5) are always aligned in the M^-plane, there 
cannot be alignments of three non vanishing vacua). The identities between kink angles hold if 

Im[l- 5(5^+^3(5,^2)] ^ Im[(53(5-^2)] ^ lm[l - 56^] _ Im[l - 5S^ -6'^ {5- (5^)] 

Re[l - 5<52 + J3(5 „ j2)] - Rc[53(5 _ ^2)] - j^c[l - 56'-] ^ Rc[l - 5^^ _ ^3(5 _ ^2)] ) 

which is satisfied if 

{S - S) (5((52 + P) -S^-P + SS {5(6^ + P + SS) - 26)) = (118) 

One irreducible component is again the abscissa axis 62 = 0. The system does not support bifurcation since all the 
minima arc now in the real axis of the complex x-plane. There are three special values: di = ±1 and di = 0; for 
61 — ±1, there exist only two minima and one kink orbit (12) and for 61 — there exist three minima, with the two 
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distinct kink orbits (12) and (23). Other possible values are: for Si < —1, there exist four minima, but only the three 
kink orbits (21), (13), and (34); for —l<Si< 0, there exist four minima, but only the three kink orbits (12), (24), 
and (41); for < (5i < 1, there exist four minima, but only the three kink orbits (14), (42), and (23); for Si > 1, there 
exist four minima, but only the three kink orbits (41), (13), and (32). 
The other algebraically irreducible component is given by 

WSl - 5Sl - - 26Sl - mlSl + WSf - 5S^ + 25S^Sl + bSlSj = (119) 

In Fig. 13 we plot the marginal stability curves which follow from the above expression, and we illustrate how the 
minima behave in the complex S plane. The two minima ±1 and any other pair of asymmetrically positioned minima 
on the curves introduces two distinct kink orbit possibilities, and so it illustrates the bifurcation phenomenon once 
again. 



1.0 



0.5 



0.0 



-0.5 



-1.0 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Figure 13: (Color online) The curves of marginal stability in the case of four minima, for e = —5. The set of (red) points 
illustrates the positions of the minima and the kink orbits for distinct values of the parameter 5 in the related regions. 

It is interesting to notice that the condition for alignment in the deformed model in the W plane, leads to (52 = 
again. But we have already seen that in this case there is no bifurcation anymore. Thus, once again we notice that 
the deformation procedure wash out the marginal stability curve, leaving no room for bifurcation also in this case. In 
Fig. 14 we illustrate the case A'^ = 4 with several distinct possibilities, driven by the complex parameter S. 

VI. FINAL COMMENTS 

In this work we have firstly dealt with the standard Wess-Zumino model driven by a complex scalar field engendering 
discrete Zn symmetry. The model is defined in terms of a superpotential, a holomorphic function of the complex 
field which contains N minima, the vacua manifold which represents a set of points with the very same Zjy symmetry 
of the model. The minima determine several topological sectors, which can be represented by algebraic curves and 
solved by first-order diff'erential equations of the BPS type, as shown in Ref. 

The main idea of this work concerns the deformation procedure developed in [Tol . [20| , which was used to deform 
the model in a way such that the set of N minima could be replicated in the entire configuration space, the plane 
described by the complex field. In this way, the algebraic orbits in field (target) space of the original Wess-Zumino 
model are also replicated in the entire plane, naturally leading to a network of defects, a spread network of kinklike 
orbits in field space. As we have seen, the idea was implemented very efficiently, and we have illustrated the procedure 
with three interesting cases, involving the Z2, Z3 and Z4 symmetries, with the sets of minima forming an equilateral 
triangle and a square in the last two cases, respectively. 

We have then investigated the more general case, with models engendering TV non symmetric minima, with = 3 
and A'^ = 4. The case iV = 3 is driven by a single complex parameter i5, and we have illustrated the results with two 
distinct values of 5, showing how it modifies the regular structure that we have obtained in the symmetric case. The 
case = 4 is more complicated, since it is controlled by two complex parameters, i5 and e. We have then made an 
interesting simplification, reducing the model to a single complex parameter, with e = —S. In this case, we have also 
illustrated the results with two distinct values of the parameter, to show how it changes the regular structure of the 
symmetric case. In the non symmetric case, we have also examined bifurcation and the marginal stability curve for 
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Figure 14: (Color online) The bifurcation curves and some related illustrations with (green) points in the first column repre- 
senting distinct values of the complex parameter 5 in the case of four minima. 



both n = 3, and ti = 4, in the last case after introducing the simphfication which leads to models described by a single 
complex parameter. In both cases, the complex parameter gives rise to a diversity of very nice possibilities, but the 
deformation washes out bifurcation from the deformed models. 

The present investigation poses several issues, one of them concerning the N = 5 and = 6 cases, which follows 
the natural course of this work. In the asymmetric case of four minima, we could also consider other relations between 
the two parameters S and e. Another issue concerns the tiling of the plane with regular polygons, which is constrained 
to appear with regular hexagons, squares and equilateral triangles, in this order of decreasing efficiency. It would 
be interesting to search for possible connection among these regular tilings, the respective basic polygons, and the 
deformation procedure. It is also interesting to consider the issue studied in Ref. (25[, in which one deals with the 
problem of marginal stability in terms of forces between soliton states in the case of three minima. A natural extension 
to the case of four minima seems desirable, and could also include investigations on how the forces between solitons 
would behave under the deformation procedure used in the present work. 

Another route of interest is related to a topological change in the field plane itself: if we let the fields to live in a 
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genus or genus 1 Ricmann surface, that is, if we consider the target space to be the two-sphere §^ or the two-torus 
= §1 X S^, instead of the plane R^, we could be able to tile the surface or with other patterns, nesting distinct 
networks of defects. This last case may possibly lead us to an issue of current interest: the conflict between geometry 
and topology, related to the geometric features of the tiling with regular polygons and the topological properties of 
the configuration space itself. It seems plausible that the modular parameter of the target space is forced to be 
our T = wa/tJi and the sphere §^ is restricted to a 2:1 embedding of this T^(ti)) in CP^. These and other issues are 
now under consideration, and we hope to report on them in the near future. 
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